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Abstract 

Several mechanisms of bar formation in stellar galactic disks are considered, including Toomre swing amplification 
and normal mode approach. On example of the well-known model of Kuzmin- Toomre using N-body simulations 
it was shown that the stellar bar is developed as a result of the evolution of an unstable normal mode. The 
pattern speed and the growth rate found agree well with linear perturbation theory. Nonlinear evolution of the 
bar is followed. Role of the growing transient spirals in bar formation is discussed. 

1 Introduction 

According to observations, bars are more common in interacting galaxies, but the percentage of 
occurrence of isolated galaxies with bar is also high. Consequently, there must be an internal 
mechanism for the bar formation. Unlike spirals, which presumably require gas inflow from 
the outside (Bournaud & Combes 2002), long-lived bars are easily obtained in the simplest 
models of stellar disks (Hohl 1971). However, despite the extensive, mainly numerical work, the 
generally accepted theory for bar formation does not still exist. To date, several mechanisms 
have been proposed, two of which will be touched here: swing amplification theory of transient 
spirals (Toomre 1981), and the unstable precessing modes by Polyachenko (2004, 2005). See also 
Lynden-Bell (1996) for some other formation mechanisms. 

One of the problems in the theory of spiral structures is the so called winding problem: 
differentially rotating disk quickly winds up and dissolves material arms. The hypothesis that a 
pattern rotates without changing its shape due to gravity has been proposed for the first time 
by Lindblad (1963, 1964), and then, in the context of stationary density waves, by Lin and Shu 
(1964, 1966). In the latter, a WKB-theory of tightly- wound spirals has been proposed. Then 
Toomre (1969) shows that such spirals evolve: propagate inward and wind up like material arms. 
So, instead of stationary density waves, the attention was switched to unstable density waves. 



In the textbook by Binney and Tremaine (2008) which is now the standard in stellar dynamics 
courses, a physical interpretation of bar-mode instability based on swing amplification mechanism 
is proposed. The mechanism is based on the effect of significant (tenfold) increase of the amplitude 
of spiral waves in the corotation region during its transformation from leading to trailing. The 
trailing waves passing through the disk center can become leading ones and re-enter the region 
of corotation. Thus multiple amplification and instability may occur. According to the theory, 
spiral patterns occur when growth rates are large, while bars occur for much lower growth rates, 
when the amplitudes of the leading and trailing waves are nearly equal. Then the interference 
pattern (in the linear regime when amplitudes are small) will have a lumpy structure with 90° 
spacing of the successive density maxima (see Binney, Tremaine 2008, Fig. 6.23). 

There are alternative modal theories in which observed structures are unstable normal modes 
of the galactic disk. In contrast to the WKB theory by Lin- Shu, where the disk is considered 
locally, here a problem of finding of the global modes that span the entire disk or its substantial 
part is posed. 

Note that calculation of the disk normal modes is a difficult task. Until recently, the only 
method was the Kalnajs matrix method (1977). Almost all available studies are adaptation of 
this method for a specific model. The method is characterized by a nonlinear equation for the 
unknown oscillation frequency, which makes it difficult to find solutions. 

Recently Polyachenko (2004, 2005) proposed a new method for the calculation of the normal 
modes, where the unknown oscillation frequencies are obtained from a linear eigenvalue problem. 
Solution of the problem are eigenvalues of a matrix, so there are no missing roots. In addition, 
this approach can be generalized to more complex systems^ and allows the insight for the bar- 
formation mechanism. 

In Polyachenko (2004) the galactic bar is an unstable normal mode in the system of precessing 
orbits. The author deduces an equation for normal modes away from resonances (except the inner 
Lindblad) and shows that stellar disks can support stable normal modes with pattern speeds 
above the maximum of precession rate, localized in the central region of the disk. These normal 
modes become unstable when the "outer" resonances (the corotation, outer Lindblad, etc.) are 
taken into account. 

Lynden-Bell and Kalnajs (1972) have shown that an outside disturbance, such as a sateUite 
galaxy, cause the ILR stars to lose angular momentum, while all outer resonances gain it^. When 

similar method for spheres was first given in Polyachenko et al. (2007); loss cone gravitational instability, 
radial orbit instability, and generalized polytropic models were studied in Polyachenko et al. (2008-2011). 

^The same is true for the energy: the change of energy and angular momentum is related by AE/AL = Clp, 
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Fi gure 1; The frequency curves Q — k/2, fi, f2 + k/2 v.s. radius R for a typical galaxy with slowly increasing 
rotation curve (solid lines, bottom to top, respectively), where fl{R) and k{R) - angular and epicyclic frequencies 
of star oscillations. Q — k/2 determines the precession rate of stars on nearly circular orbits. The positions of the 
main resonances for two values of pattern speeds Qg and fib of two-arm (m = 2) spiral pattern is shown. The 
intersection of ilg with curves indicates the position of the inner Lindblad resonances (ilLR, oILR) , the corotation 
resonance (CR) and outer Lindblad resonance (OLR). At the higher speed fib inner Lindblad resonances are 
absent. The intersections of the horizontal line fib with dashed-dotted curves show more distant resonances. 
Horizontal lines above the precession maximum show the area of "almost resonance" for fit,, see text. 

a growth rate of the perturbation is small, only stars on resonances interact with disturbance. 
However, for larger growth rates, stars in adjacent regions of the phase space are involved in 
exchange of the angular momentum. 

When the unstable mode develops in the isolated disk, a wave composed of stars in ILR region 
plays a role of the outside perturbation at outer resonances. At the same time, ILR stars are 
influenced by a perturbation of the gravitational potential from stars on the outer resonances. 
The total angular momentum of the wave always vanishes by the conservation law. 

Fig. 1 shows the location of main resonances for a system with slowly increasing rotation 
curve (for a Plummer-like model, see Section 2) for two-arm (m = 2) spiral pattern and two 
pattern speeds fls and f2b. For Qg below maximum of the precession curve Q — k/2, there are two 
inner Lindblad resonances (ILRs), the corotation resonance (CR), the outer Lindblad resonance 
(OLR), and further resonances. For fib, the ILRs arc absent. 

The observed bars in galaxies have pattern speeds fib slightly exceeding the precession curve 
maximum (Combes & Elmegreen 1993, Lynden-Bell 1996). A bar- mode with pattern speed fib 
and small growth rate cannot exist, since the only source of the angular momentum for the 
outer resonances is absent. However, for significant growth rates it is sufficient to have "almost 
resonance", from where the angular momentum can be effectively taken (in Fig. 1 the area market 
by horizontal lines), 
where flp is a pattern speed. 
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Thus, using numerical simulation, we can make a comparison of the two theories of the bar 
formation on several points: 

Process of bar formation. If the bar is formed as a result of nonlinear evolution of an unsta- 
ble normal mode, its formation can be easily tracked using Fourier analysis of the density 
distribution. On the contrary, if the bar is formed from a lumpy structure that is a su- 
perposition of transient spirals, a clear picture of the evolution of one mode should not be 
found. 

Growth rate in the linear regime. In swing amplification theory, the growth rates of bar 
modes are low (several percents of the pattern speed) . According to the theory of unstable 
precessing modes, growth rates can not be small. 

Form of the pattern. As a rule, a disk possesses several unstable modes. Evidently, the mode 
with the highest growth rate absorbs the others in the non-linear regime. Therefore, the 
theory of unstable normal modes predicts that the bar is formed from a spiral with few 
density maxima and the highest growth rate. On the contrary, swing amplification theory 
predicts lumpy structures made of superposition of the leading and trailing spirals. 

It makes sense to perform the comparison using a well-known model, such as Kuzmin-Toomre 
disk with Kalnajs (1976) distribution function (DF) rriK — 6, Lc — 0.25. This model was both 
studied numerically (Athanassoula, Sellwood 1986) and analytically (Polyachenko 2004, 2005; 
Jalali, Hunter 2005). In the latter, frequencies of unstable modes and their patterns were found 
for several models, including the mentioned above. It is worth comparing our numerical results 
with their linear perturbation analysis. Note that similar task for two-arm (m = 2), and also for 
m = 3, 4 perturbations has been made by Khoperskov et al. (2007). The connection between the 
linear bar-mode and the well- developed bar can be established using power spectra of unstable 
two-arm perturbations. 

Section 2 describes the model, gives expressions for the construction of the equilibrium model 
in presence of softening, and discusses some of the features of numerical simulations. In Section 
3, for the Kuzmin-Toomre model a comparison with the linear perturbation theory and analysis 
of power spectra for different radii and spirahty is presented. Section 4 resumes the results and 
discusses a role of the swing amplification mechanism in the bar formation. 
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2 Numerical simulation of a stellar disk 



2.1 The dynamical model 

The simplest models for study of the large-scale bar formation are self-gravitating models with 
rotation curve slowly increasing with radius. A well-known example is Kuzmin-Toomre disk with 

the surface density profile 

where a is the radial scale, is a disk mass. The corresponding gravitational potential is of 
Plummer form with the same radial scale a: 

'^^'^-^ ■ 

Numerical modeling of coUisionless systems requires gravity softening. It can be achieved, for 
example, by the replacement of the Coulomb interaction law by Plummer law: 

R + £2)1/2 ' y-') 

with a cutoff parameter e. The disk gravitational potential will be also of Plummer form: 

but with the radial scale b — a-\- e. The form of the potential is preserved if one adds a passive 
spherical component of mass with a density 



(r is a spherical radius). A full potential then requires change of to M = + Mh in 
expression (4). 

Radial density distribution (1) can be obtained with variety of distributions in phase space. 
Here we use the Kalnajs DF (Kalnajs 1976). For this we introduce new dimensionless variables 
for the energy E and potential $^: 

e^E/^^{Q) , ^ $7$^(0) = ^=L= , (6) 
and for the angular momentum L ^, and the radius R: 

x^-{-2Ef/'L/{R^^^{Q)) , y^R^^{R)/{R.^%Q))^-j^= , (7) 



^In axisymmetric system the integral is a component of the angular nionicntum along the axis of symmetry. 
In this article, we are only interested in this component. So, for the sake of brevity, we shall call it an angular 
momentum, although it is not entirely correct. 
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where = lim m'(R)/^^{0) = b. 

The surface density (1) can be represented as a function of two variables in the form 

a(w,y)^w^^r{y) , (8) 

where ttik > is a free parameter of the model, while function T{y) has a form: 

r(y)^:^,il-yT-"'-^^'(i + Sy')-'/' , S^b'/a'-l . (9) 



f{e,x)=e---'g{x) , = (27:)' J ^^^f{e,x) , (10) 



A DF in phase space for models in the form (8) can be presented as follows 

dEdL 

where L) is a frequency of radial oscillations. For zero cutoff parameter (9) coincide with 

(31) of Kalnajs (1976) and analytic expression (as a series) for the DF is possible. For £ > 0, the 
analytic expression can not be obtained and one must use the integral representation: 

g{x) = {-(2/a;)™V[2^*^(0)]} |P™K-i(l)^^K^r(a;)]- 

-^4k-i(1)^™"^(^) + /^^k-i(^)W"M^^)^4 • (11) 

J 

Primes denote derivatives of Lcgcndre functions PmK-u -PmK-i(l) = 1) ^Lk-i ~ '"^k(^k — l)/2. 

DF (10) presumes all stars to rotate in the same direction, L > 0. The surface density of the 
disk remains unchanged, if direction of rotation for some stars with angular momentum L < L(. 
is changed. A DF with retrograde stars can be written as: 



f{E, L) 



r , L<-L, , 

i(2 + 3L/L,-(L/L,)3)/ , \L\ < , (12) 



. / , L>Lc , 

where / stands for /(e, x) expressed through coordinates and velocities. 

Thus, Kalnajs models of Kuzmin-Toomre disk have three parameters (mx, q, L^): parameter 
mx characterizes the dynamic heating of the disk, q = M^/M is a fraction of the disk mass to 
the total mass in the system disk + halo, Lc regulates the amount of retrograde stars. 

2.2 The A^-body model 

Study of disk modes within the galaxy plane and subsequent nonlinear evolution of bars requires a 
two-dimensional scheme in which motion of particles is limited to galactic plane. Consequently, 
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the position of particles in the phase space is characterized by its polar coordinates R, ip and 
velocities vr and v^p. Masses of particles are considered to be identical and equal to M^/N, 
where N is a. number of particles involved in a simulation. 

The usual procedure for selecting the initial state involves random selection of energy E and 
angular momentum L for particles according to the DF 

and then random selection of radius and azimuth. Prepared in this way, the model will have a 
fairly high level of initial perturbations when number of particles N 10^. This means a short 
linear stage of the evolution insufficient for frequency analysis. More optimal way to distribute 
particles in the phase space is to choose values of E at equal increments of the mass integrated 
over L, and then values of L at equal increments of the mass along the cut at that E. Then Nr 
equally spaced particles with the same radial and azimuthal velocities are places along the circle. 

If force calculation in the azimuthal direction is limited only to the first Nr — 1 harmonics, at 
the initial stage these forces will be at the level of machine roundoff error. This scheme provides a 
quiet start for evolution of nonradial perturbations (Sellwood and Athanassoula 1986) and allows 
for a long linear growth of unstable modes. 

A PM-scheme with FFT is used for potential and force calculation. Spatial grid in polar 
coordinates (i?, </?) has n^p evenly spaced nodes in azimuth and ur nodes along the radius with 
concentration towards the center (Pfenniger and Friedli, 1993). 

For integration of stellar trajectories we use a usual leap-frog scheme. Accuracy of the initial 
distribution is checked by calculation of the virial ratio of the kinetic and potential energies, 
which should be close to —1/2 in the limit e — )■ 0. Also, the system must stay in equilibrium for 
a long time if the tangential force is set equal to zero. Accuracy of time integration is controlled 
by calculation of integrals of motion - the energy and the angular momentum. 

3 Appearance of a bar 

This section contains details of the evolution of Kalnajs model {tjik = 6, g = 1, Lc = 0.25). The 
model is well known in literature and was studied by many authors. It is characterized by a flat 
profile of Toomre parameter Q ~ 1.5 (Athanassoula and Sellwood, 1986), and moderate value of 
Julian-Toomre parameter (1966): 
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Parameters Q and X determine efficiency of swing ampfification mechanism (Toomre 1981). The 
Plummer gravitational potential gives the slow growth of the rotation curve in the center, so 
the ILR doesn't interfere the passage of the wave packets. Therefore, this model is ideal for the 
emergence of swing-amplified unstable modes. 

3.1 The normal bar-mode 

Fig. 2 shows frames of evolution of the model with cutoff parameter e — 0.2. It is assumed that 
a = M = G = 1, that gives for time units ~ 10 {a/Wi^cf/'^iM/W^MQ)-^/'^ Myr. The grid 
was riR — — 64, step of time integration h — 0.024533, Fourier expansion of azimuthal forces 
was cut at Nh — 9, number of particles on circles Nr — 10. Accuracy of energy and angular 
momentum conservation was 6e ~ 10~^ and Sl ~ 10~^. Full time of the evolution Tmax = 245.33. 

The evolution of the stellar disk can be divided visibly into three stages. In the first stage 
t < 120 the image of the disk remains practically unchanged preserving an axisymmetric form, as 
in frames t = 75, 100. This is followed by a relatively short second stage of formation of a pattern 
(bar-mode). The third stage t > 150 is characterized by strongly non-linear bar evolution, when 
bar changes slowly in length and its rate of rotation. 



Time = 75 Time = 100 Time = 125 




-3 -2 -1 1 2 3 -3 -2 -1 1 2 3 -3 -2 -1 1 2 3 



Figure 2: Frames of the evolution of the Kalnajs (6, 1, 0.25) model with cutoff parameter e = 0.2, for moments 
t = 75, 100, 125, 150, 175,245. Number of particles N = 10^. For a = Ikpc and M = lO^M© the last frame 
T = 245 corresponds to ~ 2.45 Gyr. 
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Figure 3: The dependences of ln|A2(p, t)| v.s. spirality p shown for different moments of evolution t = 
0, 10, 20, 240. Crosses show maxima for each curve. 

A more detailed view of what is happening in the disk at these stages can be obtained by 
Fourier analysis and the expansion in logarithmic spirals (Sellwood and Athanassoula 1986). 
Consider the expansion coefficients 

Am{p,t) = J J S(r, t) exp(— i[my9 + plnr])rdrdv9 . (15) 



The pitch angle of log spirals is constant and equals arctg(m/j9), where p is spirality (trailing 
spirals have positive spirality, bars correspond to p = 0). In what follows, we are interested only 
in the two-arm perturbations, i.e. m = 2. Having the A^-body model surface density in the form 

S(r,y,,t) = ^J]5[r-r,(t)]5b-^,(t)] , (16) 

i 

where rj(t) n <fi{t) are coordinates of particles, one obtains 

A„(p,t) = ^exp(-z[mv9i +plnri]) . (17) 

i 

Fig. 3 shows the curves (j9, t) \ characterized by rather broad plateau, indicating the presence 
of components with different p in the perturbation. At the beginning the plateau is localized 
in the interval — 1 < p < 3 with average spirality p = 1. Then, in the course of the pattern 
development, it moves left in the interval —2<p<2 with the average value p = 0. Crosses, 
indicating the positions of the maxima of the curves also show a tendency of the pitch angle 
increase with growth of the disturbance amplitude. It is important, however, that function 
|A2(p, t)| has no pronounced maxima at large and opposite values of p corresponding to leading 
and trailing tightly-wound spirals, as can be expected in the Toomre theory of bar instability. 
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Fig. 4 shows time dependence for T{t) = max In |A2(p, t)| and complex phase A2{p,t) = 

p 

\A2\e'^'^ found at fixed t for values of p for which In |A2(p, t) \ peaks. The smoothness of the curves 
and the presence of broad plateau in Fig. 3 indicates that the result should not depend on a 
particular form of spirals used for decomposition. 

The main thing to note in Fig. 4a that curve r(t) follows precisely the linear law at t < 140, 
up to amplitudes In 1 — —2.3. The linearity means that the amplitude of the perturbation 
grows exponentially with time, i.e. oc exp{'yt), where 7 is the growth rate. The slope of the linear 
region of r(t) determines the growth rate of the instability. After reaching a level of about the 
unperturbed surface density at t ~ 170, the bar amplitude saturates at the level In IA2I ~ —1. 

The phase (j){t) demonstrates even better linear behavior at t < 150, see Fig. 4b, which means 
a constant rate of rotation of the global pattern 

«5 = -f . (18) 

^ m at 

After the linear stage, the rate of rotation slows down. 

A fairy long period of linear evolution t < 140 with constant rotation velocity and growth 
rate suggests the existence of an unstable normal mode. Note that no sign of the pattern can 
be visibly detected for most of this period (up to about t < 117, In |y42| < —4), and a quite well 
shaped bar appears only at the end of the linear phase at t ^ 150. 

The pattern speeds and growth rates calculated for different values of the cutoff parameter e 
are shown in Fig. 5 and in the table. The least square quadratic extrapolation to £ = gives for 
the pattern speed without softening ^ 0.428, and for the growth rate 7*^ ~ 0.312. 

Jalali and Hunter (2005) independently investigated the model without softening by linear 
perturbation analysis. They have found two unstable normal modes. For the most unstable 




Time Time 

Figure 4: Time dependences of r(t) = niacin \A2{p,t)\ (a) and complex phase 0, A2{p,t) = \A2\e^'^ (b) found 

p 

at fixed t for values of p for which ln|v42(p, i)| peaks. Dashed lines show linear behavior for both dependences 
obtained by least square method for 10 < t < 130. 
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£2 (0) =0.428 



C! 0-35 - 



b) 



Y(0) =0.312 



Figure 5: Pattern speeds fip and the growth rates of the mode 7 at values of cutoff parameter e = 0.2, 0.15, 
0.1, 0.07, 0.05 and its extrapolation to e = 0. 

mode, the pattern speed Vtp ~ 0.445, and the growth rate 7 ~ 0.308, for the second unstable 
mode - f2p ~ 0.294, and the growth rate of 7 ^ 0.109. Our values are very close to the values of 
the first mode found by Jalali and Hunter. This once again confirms the conclusion that we are 
dealing with the most unstable normal mode. 

A shape of the mode remains unchanged in the linear stage of the evolution. It looks like an 
open spiral with one maximum (see Fig. 3a in Jalali, Hunter 2005). In the next section we will 
see how this unstable bar mode is transformed into a bar in the nonlinear stage. 

3.2 Nonlinear evolution of the bar mode 

The numerical experiment demonstrates the presence of one unstable mode. The developed bar 
on the late stage of the evolution is a nonlinear structure with an amplitude comparable with 
the axisymmetric background. Examination of the power spectra of disk oscillations at different 
time intervals allows one to associate a small amplitude unstable bar-mode and the non-linear 
bar. 

Let's introduce two kinds of power spectra for different spirality p and radii in the interval 
+ s). First of all, for determination of pattern speed Vtp (or real part of the frequency 
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0.20 


0.15 


0.10 


0.07 


0.05 
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0.229 


0.267 


0.314 


0.344 


0.366 


0.428 


0.445 


7 


0.076 


0.110 


0.161 


0.196 


0.229 


0.312 


0.308 



Table 1: Pattern speeds Vtp and the growth rates 7 of the mode found from slopes of T{t) and 
0(t) at values of cutoff parameter e = 0.2, 0.15, 0.1, 0.07, 0.05 and its extrapolation to £ = 0. 
The last column contains figures of independent linear perturbation analysis by Jalali and Hunter 
(2005) for the model without softening. 
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Reu; = mflp) it is necessary to exclude the growth of disturbances. One way of achieving it is to 
exclude the exponential growth factor exp(7f ) from the amplitude A2{p, t'). Then the expression 
for the spectra will be limited to the linear stage, while the growth rate is constant and equal to 
7 = Imco. A more general way is to normalize the amplitude to the maximum of its modulus 
r(t'), thus the expression for the spectra can be used both in linear and nonlinear regimes, even 
when there is no increase in amplitude. Therefore we define the power spectrum, corresponding 
to different values of spirality, by the expression: 

2 



t+s 



Mp,t' ) , 

m 



e'"*"^* Hl{t')dt' 



(19) 



where H^{t') is the Hann filter: 
Similarly, for the radial spectrum: 



1 — cos 27r 



t' -t 



(20) 



t+s _ 



^^^^^P e^rn^vt' Hf(t')dt' 



(21) 



where 



A2{R,t) — j J Ti{r,(f,t) exp{—imip)rdrd(p — exp{—i'mipi) . (22) 

R-AR * 

Prime in the last expression means summation over particles in the ring {R — AR, R + AR). 

Fig. 6 shows power spectra for the model with cutoff parameter e — 0.2 for hnear {25 < t < 95) 
and nonhnear (175 < t < 245) phases of bar evolution. 

The peaks of the power spectra show the position of the coherent structures, in which a large 
number of particles are involved. Fig. a) and b) are the power spectra corresponding to different 
values of flp and spirality p, calculated from (19). Similar to frames of the evolution displayed 
in Fig. 2, one can trace time evolution of the power spectra. During the linear phase, contours 
of power spectra keep their position near the point ~ 0.23 and spirality p — 1, corresponding 
to an open spiral with an average pitch angle a — 60°. During the transition period and the 
nonlinear stage, one can observe the pattern's slowdown and its transformation into a bar (with 
p^O). 

Figures b) and d) show maxima of the radial power spectrum (21) revealing radial localization 
of the mode. The frequency curves taken from Fig. 1, but with softening, allow one to compare 
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Figure 6: Isolines of the power spectra (19) (figures a,c) and (21) (figures b,d), for the linear phase of the 
evolution of the bar-mode 25 < t < 95, and for the nonlinear evolution of the bar 175 < t < 245. Filter width in 
all cases s = 70. Outer contours on each figure correspond to half of the spectra maxima, next contours increment 
5% of the maxima. Solid lines in figures b), d) show frequency curves il — k/2, fl, ft + k/2, similar to Fig. 1. 

the position of the spectrum maxima with the position of the maximum of the precession curve 
Q — k/2. It can be seen that the disturbance is concentrated over the maximum, according to 
the theory of bar formation by Polyachenko (2004). After the slowdown, the bar speed becomes 
shghtly higher than the maximum of the precession curve. 

The described temporal behavior of the pattern speed Qp and spirality p is clearly seen in 
Fig. 7, that shows location of maxima of Vf{p, Qp) versus time. The plot in Fig. 7a is consistent 
with our findings based on analysis of Fig. 4b, about the behavior of the pattern speed, taking 
into account a time lag s/2 of the filter (20). 

4 Discussion 

The quiet start technique by Sellwood and Athanassoula (1986) allows one to construct stellar 
disk models with low level of initial perturbations and to study their linear and nonlinear evolution 
using different Fourier spectra. 

Our investigation of Kuzmin-Toomre disk model with flat Toomre parameter proflle Q ~ 
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Figure 7: Time dependence of the pattern speed Up{t) (a) and spirality p{t) (b), found by location of power 
spectra maxima (19) and (21), for the bar-mode (cutoff parameter e = 0.2). The filter width is keep constant, 
s = 30. 

1.5 reveals the presence of a disk unstable normal mode. This mode dominates in the disk 
from the beginning of the simulations. A comparison with the independent results of the linear 
perturbation analysis by Jalali and Hunter (2005) shows that it is the most unstable normal 
mode of the stellar disk, the surface density profile of which has a unique maximum. 

Upon reaching the level comparable with the axisymmetric background, the mode is governed 
by non-linear evolution: the amplitude stops growing and the pattern speed slows down. Using 
power spectra we see that the pattern preserves the open spiral shape during the linear evolution, 
then it transforms into a more open spiral and eventually a bar. 

Frequency analysis of perturbations hasn't shown any other unstable modes. Clearly, this is 
the result of our crude analysis that removes all the modes except the most unstable one. 

The model meets necessary conditions for effective swing amplification of transient waves. 
Thus, it would be natural to expect the presence of leading and trailing spirals that form the 
bar. However, we could not find it neither visually nor using a more detailed frequency analy- 
sis. According to Binney and Tremaine (2008), in case of considerable growth rates the swing 
mechanism should give spiral patterns, while bar-modes are characterized by low growth rates 
and lumpy structures. We have seen that it contradicts to results of our numerical experiments, 
where the bar is formed from the most unstable mode with one maximum. 

Note that in the literature there are different interpretations of bar instability even in the 
framework of the swing amplification. For example, in the recent review by Sellwood (2010) 
one can find: «Toomre (1981) provided the most important step forward by elucidating the 
mechanism of the bar instability (see also BT08 sect. 6.3). Linear bar-forming modes are standing 
waves in a cavity, akin to the familiar modes of organ pipes and guitar strings. Reflections in 
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galaxies take place at the center and at the corotation radius, except that outgoing leading spiral 
waves incident on the corotation circle are super-reflected into amplifled ingoing trailing waves 
(i.e. swing amplification), while also exciting an outgoing transmitted trailing wave. The feedback 
loop is closed by the ingoing trailing wave reflecting off the disk center into a leading wave, which 
propagates outwards because the group velocity of leading waves has the opposite sign to that 
of the corresponding trailing wave. The amplitude of the continuous wave-train at any point in 
the loop rises exponentially, because the circuit includes positive feedback.* 

First, the reflection takes place from a Q-barrier located near the corotation. A wave cannot 
propagate freely under the barrier. Second, the outgoing transmitted trailing wave behind the 
Q-barrier hasn't been mentioned in Toomre (1981) at all. For the first time it appeared in waser- 
I mechanism by Mark (1974) for spiral disturbances. Meanwhile Sellwood (2010) repeats the 
so-called waser-II mechanism for bar- modes proposed by Bertin (1983), which is the example of 
a modal approach, where bars are treated as unstable normal modes. 

Waves outside corotation {Q < Qp) have positive angular momentum, while waves inside 
corotation {Q > Qp) have negative angular momentum. Penetration through the Q-barrier 
provides a mechanism for the instability since generation of positive angular momentum wave 
strengthens a negative angular momentum wave and vise versa. The presence of a third wave 
is important, since it explains the source of instability by means of conservation laws. If swing 
amplification mechanism does not require the third wave, it is quite different from what is usually 
meant. 

It seems, however, that Toomre's approach to bar formation in galaxies involving transient 
spirals and due to swing amplification is not a necessity. There are many reasons, especially the 
relatively recent observations in the near-infrared, to believe that instead barred spiral structures 
should be naturally interpreted in terms of global modes. The amplification mechanism can 
be the same independently whether the structures are transient or not. The transient-swing 
amplification picture is a mathematical description that has a direct counterpart in the dynamics 
of steady wavetrains at corotation. 

Comparison of our theory of bar formation with waser-II mechanism requires further investi- 
gation, since both mechanisms are modal. For example, one can investigate the effect of feedback 
break on the stabihty properties of stellar disks and what is the relative role of the corotation 
region, the outer Lindblad resonance, and other external resonances in bar mode instability . 
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